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ABSTRACT 

Using the recently constructed basis for local operators in free SO(N) gauge theory we derive 
an exact formula for the correlation functions of multi trace operators. This formula is used 
to obtain a simpler form and a simple product rule for the operators in the SO(N) basis. 
The coefficients of the product rule are the Littlewood- Richardson numbers which determine 
the corresponding product rule in free U(N) gauge theory. SO(N) gauge theory is dual to a 
non-oriented string theory on the AdSs x HP 5 geometry. To explore the physics of this string 
theory we consider the limit of the gauge theory that, for the U(N) gauge theory, is dual to 
the pp-wave limit of AdS§ x S 5 . Non-planar unoriented ribbon diagrams do not survive this 
limit. We give arguments that the number of operators in our basis matches counting using 
the exact free field partition function of free SO(N) gauge theory. We connect the basis we 
have constructed to free fermions, which has a natural interpretation in terms of a class of 
|-BPS bubbling geometries, which arise as orientifolds of type IIB string theory. Finally, 
we obtain a complete generalization of these results to Sp(N) gauge theory by proving 
that the finite iV physics of SO(N) and Sp(N) gauge theory are related by exchanging 
symmetrizations and antisymmetrizations and replacing N by — N. 
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1 Introduction 

In a previous paper [I] we have initiated the study of local operators in SO(N) gauge theory 
which have a bare dimension that can depend parametrically on N. For these operators, one 
needs to sum more than just the planar diagrams to capture the large iV limit. We dealt 
with this problem by employing group representation theory to define local operators which 
generalize the Schur polynomials of the theory with gauge group U(N). We found that the 
free field two point function is diagonalized by our operators. In this article we will extend 
our understanding in a number of important ways. 

A basic result of [T] is the basis of local operators, given by 

Or(Z) = T^y E XR^a^^t^Z^Z^ ■ ■ ■ Z*-*» (1.1) 



1 



For SO(N) with N even, we also need to include^ 



Qr(Z) 



(iV + 4)! 



Jli2---«JV-l2jV*JV+l*iV+2«JV+2*iV+l'"«JV+p-l«JV+p*JV+p*JV+p-l 



Jl]2---jN+2p-ljN+2p 
X 73lh 733H . . . ^3N+2p-i3N+2p-2 %3N+2p-l3N+2p 
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We will focus our discussion on the Or(Z) which we understand better than the Qr(Z). 
The operator label R for Or is a Young diagram with 2n boxes, that is, an irreducible 
representation of the symmetric group S^n- To obtain a non-zero operator Or(Z), n must 
be even. Thus, 2n is divisible by 4. In fact, the only representations R which lead to a 
non-zero operator are built from the "basic block" Lu as we now explain. Choose a partition 



of n/2, or equivalently a Young diagram r with n/2 boxes. The set of valid representations 
R are obtained by replacing each box in r by the "basic block" Hi . To reflect the relation 
between r and R, we use the notation r = R/4. Thus, the number of gauge invariant 
operators of the type Or(Z) built using n fields is equal to the number of partitions of |. 
As an example, for n = 4 the allowed labels are 



# 2 /4 
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Ri 



Ri 
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Given the form of the Young diagrams R we consider, we will list the row lengths of R as 
(r 1; r 1; r 2 , r 2 , • ■ ■ , rjv , rw). In section 2 we start by using the results of pQ to derive an exact 
formula for the free field theory correlation functions of multi trace operators. The result is 
given in equation (12.1 5p . Using this formula we explain in section 3, how to obtain a simpler 
form for the Or(Z) of our SO(N) basis. The resulting operators, with a more convenient 
normalization, are 



X°s(Z) 



(f)«; 



eSr 
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With the new normalization, the two point function of these operators is 



(xKzh^z)) = 5 RS J] r, 

jgodd boxes in S 



1 Recall that there are two invariant tensors for SO{N): the Kronecker delta and the e^i 
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We can 



use either of these tensors when contracting indices to obtain gauge invariant operators. For SO(N) with N 



odd e ili2 ... il 
a single t ixi 



has an odd number of indices, so that we can't build a gauge invariant operator that uses only 
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A few comments are in order. Each box in a Young diagram can be assigned a factor, denoted 
Ci in the above formula. A box appearing in column a and row b has factor N + a — b. The 
right hand side is equal to the product of the factors of the boxes in every second row as 
shown below 
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We called these the odd boxes in pQ because they referred to boxes labeled with an odd 
integer in a Young- Yamanouchi labeling of the states in the 5*2 n irreducible representation. 
With the new simplified form of the SO(N) operators, we are able, in section 3, to give a 
product rule for our operators. The product rule is (S/A h R/4 h ^) 

X°s(Z)Xr(Z) = 9R/±s/4t/4Xt(Z) (1-8) 

T/4h ni+n2 

where gB,/4S/AT/i is the Littlewood- Richardson coefficient. 

These results constitute a rather complete understanding of the local operators in SO(N) 
gauge theory, comparable to what has been achieved for the U(N) theory. This program 
was initiated in the context of U(N) gauge theory, by Corley, Jevicki and Ramgoolam in 
[2]. In particular, [2] showed that the half-BPS operators constructed using a single complex 
matrix can be described using Schur polynomials and they demonstrated that the Schur 
polynomials diagonalize the free field two point function. The study of the finite N physics 
of U(N) gauge theories is by now well developed. There are a number of bases of local 
operators that diagonalize the free field two point function [31 IU El El CO El El EE] and we 
know how to diagonalize the one- loop dilatation operator [TJ] H2] for certain operators dual 
to giant graviton branes[T3| [HI [151 HS1 HZ]- This diagonalization has provided new integrable 
sectors, with the spectrum of the dilatation operator reducing to that of decoupled harmonic 
oscillators which describe the excitations of the system [T4 T [TB I 116] . Integrability in the planar 
limit was discovered in [TS1 [JJJ] and is reviewed in [20J. For a study of the SU(N) theory see 



Given the results we have developed, we are now in a position to probe the finite N 
physics of SO(N) gauge theory. Recall that according to the AdS/CFT duality [221 E21 
El], finite N physics of the gauge theory [25] corresponds to non-perturbative (in the string 
coupling) physics of objects such as giant graviton branes [26 | 127 1 128] and the stringy exclusion 
principle [29J. M = 4 super Yang-Mills with SO(N) or Sp(N) gauge group is dual to the 
AdS 5 x ftP 5 geometry [3D]. In this case one expects a non-oriented string theory so that 
the study of non-perturbative stringy physics, which is captured by the finite N physics 
of the gauge theory, is likely to provide new insights extending what can be learned from 
the AdSsxS 5 example which involves oriented string. For studies in this direction see [3Tj . 
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Computations in the gauge theory, must sum both the planar and the non-planar diagrams. 
At the non-planar level there are genuine differences between the U(N) and the SO(N) or 
Sp(N) gauge theories. Recall that matrix model Feynman diagrams in double line notation 
represent discrete triangulations of Riemann surfaces [32J. For Hermitian matrices we deal 
with oriented triangulations whereas for symplectic or anti-symmetric matrices unoriented 
triangulations [55]. In general, a Feynman diagram is weighted by 

y2g-2+b+c N -c-2g+2^ /-^gX 

where N is the number of colors, A = gy M iV is the 't Hooft coupling, g is the number 
of handles, b the number of boundaries and c the number of cross-caps on the surface. 
Thus, for the SO(N) or Sp(N) gauge theories, the leading non-planar corrections come from 
ribbon graphs that triangulate non-orientable Feynman diagrams with a single cross-cap. 
The large N limit of correlation functions of operators with a bare dimension that depends 
parametrically on N are sensitive to this non-planar structure of the theory. With the goal 
of probing this structure, in section 4 we use our technology to compute free field theory 
correlation functions of multi trace operators, to all orders in 1/N. Given these correlators, 
we can consider the double scaling limit defined by [3H [35] 

N — > oo and J — > oo with — fixed, gyu fixed (1-10) 

where J is the number of fields in the gauge theory operator. In this limit some non- 
planar diagrams (string interactions) survive, giving a non-trivial normalization of the two 
and three-point correlators of single trace operators. This limit is particularly interesting 
because in the dual gravity picture it corresponds to taking a pp-wave limit of AdS^ x S 5 , 
a background in which the superstring theory can be quantized. We find that non-planar 
unoriented diagrams in SO(N) gauge theory do not survive this limit. 

We have argued that Or (and for N even, the Qr) give a basis. A weak point in our 
argument is that we have not demonstrated that these operators are a complete set. This 
issue is considered in section 5. We focus on N even, which is the more involved case, as a 
consequence of the fact that we may use e^^...^ when constructing gauge invariant operators. 
By counting the number of operators we have constructed, we are able to reproduce the exact 
free field partition function of the SO(N) gauge theory |36j for iV = 4, 6 and, using our results, 
we give a conjecture for the free field partition function at any even N, in (I5.3P and at any 
odd N in (jOijl . 

The basis that we have constructed allows us to study the dynamics of the gauge in- 
variant observables of a single matrix model. Of course, this problem can be reduced to 
eigenvalue dynamics which is itself equivalent to the dynamics of free fermions in an exter- 
nal potential [37]. The Schur polynomial basis for the U(N) gauge theory has a very direct 
link to free fermion dynamics. It is natural to ask if there is a similar connection between 
free fermions and the basis we have constructed. We develop this link in some detail in 
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section 6 and show that there is indeed a natural connection to free fermions. Our operators 
can be mapped to states of fermions moving in a harmonic oscillator potential, with definite 
parity and maximum angular momentum for a given energy. 

Although we will not do so in this article, we have developed enough technology that it 
would be natural to initiate a systematic study of the dilatation operator in non-planar large 
N limits of the SO(N) and Sp(N) gauge theories. A detailed study of the planar spectral 
problem of Af = 4 super Yang-Mills with gauge groups SO(N) and Sp(N) has been carried 
out in [38]. The essential difference between the theories with gauge groups U(N) or SO(N) 
is that in the SO(N) case certain states are projected out. Thus, the planar spectral problem 
of the SO(N) theory can again be mapped to an integrable spin chain [38]. It is interesting 
to ask if the new integrable sectors discovered in [HI [151 EE] are also present in large N but 
non-planar limits of SO(N) and Sp(N) gauge theory. 

It is well known that there is a close relationship between SO(N) group theory and Sp(N) 
group theory. These relations imply that the dimension of a given irreducible representation 
of Sp(N) is equal to that representation of SO(N) with symmetrizations exchanged with 
antisymmetrizations (i.e. transpose the Young tableau) and N replaced by — iV[39] . The 
QCD loop equations for SO(N) gauge theory and Sp(N) gauge theory in 3 + 1 dimensions 
enjoy the same connection [40J. This same relation has been observed in two dimensional 
Yang-Mills theory [4"T]. Motivated by this background, we proceed to study the finite N 
physics of Sp(N) gauge theory. We are able to argue that precisely the same connection 
relates the finite N physics for the orthogonal and symplectic gauge theories. In this way, in 
section 7 we obtain a rather complete description of the finite N physics of the Sp(N) gauge 
theory. 

In section 8 we outline some open problems that we find interesting. 

2 Multi trace Correlators 

Our goal in this section is to give an exact formula for the free field theory correlation 
functions of multi trace operators in SO(N) gauge theory. In our discussion below n, | and 
2n will enter at various points. The reader is encouraged to keep in mind that our operators 
are built using n fields. 

We start by choosing any partition v h ^. The partition is then translated into the cycle 
structure of a permutation. For example, if ~ = 4 there are 5 possible choices for u, namely 
(4), (2) 2 , (2) (l) 2 , (3) (1), or (l) 4 . We use these partitions v (see [lj) to construct a iu 6 S 2n 
with cycle structure given by multiplying each of the parts of v by 4. For the above list of 
partitions v the cycle structures of the corresponding permutations a iu are (16), (8) 2 , (8)(4) 2 , 
(12) (4) and (4) 4 . We will now explain, by providing a few examples, how we associate to 
each cycle structure a canonical permutation. The cycle structure (12) (4) is associated to 
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the cycle 

a 4u = (1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12)(13, 14, 15, 16) (2.1) 

while (8) (4) 2 is associated to 

a Au = (1, 2, 3, 4, 5, 6, 7, 8)(9, 10, 11, 12)(13, 14, 15, 16) (2.2) 

So, the rule for obtaining the canonical permutation is to populate the largest cycles first, 
starting from 1 and counting up, until the permutation is completely determined. Notice that 
the canonical permutations are composed of cycles with cycle lengths that are a multiple 
of 4 and further, they always take an even number to an odd number. These canonical 
permutations can be used to define a "contractor" as follows 



J2p3a(2p) 
P =l 

n 

P =i 

To see that all indices appear on the right hand side and no index appears more than once, 
it is useful to remember that o"4„ always takes an even number to an odd number. Using the 
contractors, we can define the operators 

1 /3es 2n 

These operators are particularly convenient for the question of correlation functions of multi 
trace operators. Indeed, using the identity 

= 72^1 E d *x*W ^ 

we easily find 

(^,(^)£ w -"^5* = E d *°s to m 



R 



Thus, any multi trace operator can easily be written as a linear combination of the R 
Introduce the notation 



TW(K^T f ) = {Z 2 t m (Z 2 t w •••( Z2 )l%) (2.7) 

Every multi trace operator can be written in this way for a suitable v. Thus, to evaluate an 
arbitrary two point function of multi trace operators, all we need to do is to compute the 
correlator (O^ 4 " '0^ 4m ) . Using the results in [lj, it is straight forward to see that 

x r)l9 ra 

(ot(z)o^(z)) = j2 ^{P[AF R mcrmci» (2.8) 



We can define a permutation a UfJ- by the condition 



(2.9) 



By considering an example at this point, it will be clear that a Ufl defines a unique element 
of the double coset H± v \ Szn/H^ where (H^) is a stabilizer of (a^) respectively. 
Indeed, for 



(74, = (1,2, 3, 4, 5, 6, 7, 8) (9, 10, 11, 12) 



(2.10) 



and 



a 4M = (l,2,3,4)(5,6,7,8)(9,10,ll,12) 



(2.11) 



we have 



a v „ = (1,5) 



(2.12) 



See Figured] for an illustration of this example. 




/R\ 



}4n 



Figure 1: Relating Cf and Cf with a permutation. 



With the use of a upL we can now write 



(2n)\d 
1 >' R ipes 2n 



(2.13) 



Since we are contracting row and column labels of ip in the same way, the computation 
becomes very similar to the computation of [I]. In the next formula we introduce an S n [S2\ 
subgroup that belongs to the stabilizer of Cj M . The reader should recall[T] that P[a] is defined 
using the embedding of S^^] that stabilizes (1, 2)(3, 4) • • • {2n — 1, 2n). It is now straight 
forward to obtain 



(O^(Z)O^(Z)) 



E E Tr(P M r«(^^ 2 ))^(^i^)M, 
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n\2 



(2n)!d 



Vies. 



(2n)!d, ^^WW n c * 

iSodd boxes in R, 



= §Rs2 ^)+mM^^ JJ c . (2.14) 

D 

21 i£odd boxes in R 

which is a remarkably simple formula. To obtain the last line we have used the results of 
[12] as explained in Appendix C of [TJ. Using this in a completely straight forward way we 
find 

(Tr(M^T f )Tr(K^T f )> = E ^^Xr/aMxr/M \{ Q (2.15) 

ij/4hf ieodd boxes in R 



3 A simpler description of the SO(N) basis 

In this section we will argue that the result (I2.15P allows us to write a simpler description of 
our basis, that is closely related to the Schur polynomial basis of the U(N) theory. Towards 
this end, we begin by computing the correlator (OnOg 4 "). First, it already follows from the 
results of [TJ that 

(0 R 0^) oc 5 rs (3.1) 

so that we only need to compute (OrO r 4v ). Now, noting that Or = O a R iv with v = it, we 
can apply (12.141) to find 

(0 K 0-) = Wf^)^ /4 (,i) J] c % (3.2) 

° igodd boxes in _R 

Now, since the Or constitute a basis and since (I3.2p gives the two point function of CT s il/ 
with any Or, it is clear that (13. 2p can be used to determine O a s iv as a linear combination of 
the Or. We find 

Ot = ~, Xr/^)Or (3.3) 

a R/A 

Thus, we can now write 

Tt v ® n {v{Z 2 )^) = E <W 

= E ^ 3 - 4 ) 



R/4h 



n ^«/4 



Notice that, using character orthogonality, we can now invert this relation. Indeed 



2- 1 ^ Xs/ ,(u)Tt v ^(u(Z 2 



di 



ueSr 



R/Ah 



Y J^ 2 ' Y XS/4{v)XR/4 



.V) 



ueSr 



Oi 



Y ^-f[^)!,W 

/?/4hf ^ 



ds/4 



\2 
O s 



O, 



Consequently 



O s {Z) 



^f 2f lk Y 2-^/ 4 (-)Trv-(^ 2 



(3.5) 



(3.6) 



The normalization in f !3.6p looks rather unnatural. From now on we will adopt a new 
normalization, given by 

1 



£ 2 ~ 



«(!/) 



XS/4 



z/)Tr 



(3.7) 



Notice that we continue to label our operators by 5, not by 5/4. This deserves a few com- 
ments. In the context of SO(N) gauge theory, the Wick contractions are a sum over elements 
of S 2n - To construct our operators, we have constructed projectorspQ using representations 
S of S2n- Orthogonality of our operators then follows as a consequence of the fact that these 
projectors commute with the Wick contractions and are mutually orthogonal. The two point 
function is given in terms of a product of factors of boxes in S, which were obtained [I] by 
evaluating the action of Jucys-Murphy elements on states in the carrier space of S. Clearly, 
S summarizes information about the group theory used to construct our operators: it is the 
representation that organizes the 2n indices of the Z %: > fields that appear in Os- Thus for 
example, cut offs due to the stringy exclusion principle [291 EH] cut S off at iV rows. 5/4 on 
the other hand, plays the role of a useful auxiliary label. Indeed, 5/4 is particularly useful in 
enumerating the possible operators (all possible 5/4 are allowed; only special representations 
5 are allowed) in our basis. Note also that for the operators Qr(Z), since the rows in R all 
have an odd length, there is no notion of the R/4 label as defined above. With this new 
normalization we obtain the lovely answer 



(Xr(Z)xs(Z)) 



(3. 



^ RS II Ci = ^RsfR,od<i 

tgodd boxes in S 

The formulas of this section indicate a very interesting interplay between Young diagrams 
5h2nand5/4hf. 

One immediate application of the new formula (13. 7p is in the derivation of a product rule. 
Indeed, for 5 h 2m and R h 2n we have 

1 



X°s(Z)x°r(Z) 



\ Yl Yl 2 



Xs/4(^)xW/i)Tr y «n(z/(Z T 2 )TT V ®n(ij,(Z 
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IT U — E — 



Tr y8 n +m (a(Z 2 )®f) 

/n V /my fn+m TT E E E £ ^ XS/A^XR/MXT^ O 

\2j-\2j-\ 2 J' veSm fi£Sn <T<=S n + m T/4h n + m 

/] 9r/as/4T/aXt( z ) ( 3 

T / 4 |_n±m 



where we have used the formula 



E Z) xs/4y)xR/MxT{ii 1 ° v *) (3.10) 



V 2/ ' V 2 / ' ,/eSm uGSn 



for the Littlewood- Richardson coefficient. 



4 Multi trace Correlators Again 

The answer (I2.15P gives a complete description of correlators in the trace basis. To obtain 
explicit answers we need to evaluate symmetric group characters. We explain how this 
evaluation is carried out in this section. 

In the U(N) case, very similar formulas have been obtained, exploiting the relation 
between operators written in the trace basis and the Schur polynomials. We could also 
follow this route given the new form of our operators in (13. 7p . Following this route, we 
would write multi-point correlators in terms of products of our operators and then evaluate 
these products using the Littlewood-Richardson coefficients. This approach computes the 
general multi-trace operators knowing nothing more than a character for an n cycle of S n in a 
Young diagram labeled by a hook representation and the Littlewood-Richardson coefficients. 

For SO(N) we can employ the formula (I2.15p which requires the computation of char- 
acters, beyond the character for an n cycle of S n in a Young diagram labeled by a hook 
representation. As we explain below, the evaluation of these characters is straight forward. 
Given the values of the characters we obtain below, the current computation could also be 
used to give an alternative derivation of the known U(N) correlation functions. 

To start, consider the computation for correlators of the form 

<Tr(Z^)Tr(^)Tr(^)> = 2 3 £ Xr/MXr/M U c t (4.1) 

R/AhJs i£odd boxes in R 

We have J3 = J\ + J2 and \i is a J3 cycle while v is the product of a J\ cycle and a J2 cycle. 
We know that Xr/aIh) wm om Y be non-zero when R/A is a hook. Thus, we will consider 
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only Young diagrams R/A of the formal [k, l Js k ]. A formula that will be useful i 

(N + 2k-2)\ 
(N - 2J 3 + 2k - 2)! 



n c= ( N + 2k ~ 2 V (42) 



i Sodd boxes in R with R/4=[k,l J 3~ k ] 

We know that 

X[M J a-*](M) = (-l) J3 ~ fe (4-3) 

Now consider X^i^-^^)- According to the Murnaghan-Nakayama rule[l3], this character 
is equal to a sum over all ways of extracting a border strip tableau of length J\. For each 
possible extraction we have to multiply by (— l) h where h is the height (= number of rows) 
of the removed border strip multiplied by the character of a J 2 cycle in the irreducible 
representation labeled by the Young diagram obtained by removing the border strip from 
[k, l j3_fc ]. Thus, (in the following (J 2 ) denotes a J 2 cycle) 

X[k,i-h-k](v) = -9(k > Ji)xik-Jx,i J 3- k ](( J 2)) + (-l) Jl 0(J 3 -k> Ji)xik,i J 3- k - J i](( J 2)) 

= -9(k > Ji)X[*-j 1 ,i'3-*]((^)) + (-l) Jl 0(J2 > k)x [k ,i^]((J2)) 

= 9(k > Ji)(-l) j3 " fc+1 + 9(J 2 > k)(-l) j2+Jl - k 

= 9{k > Jx)(-l) j3 - fc+1 + 9(J 2 > k){-l) Js - k (4.4) 



where 



9(k > Ji) = 1 if k > Ji 

= otherwise (4.5) 



9(J 3 -k>Ji) = 1 if J 3 - k > Ji 

= otherwise (4-6) 

It is now clear that 

X[fe,i J 3- fc ](/^)X[fe,i J 3- fe ]( z/ ) = 9(J2 >k)- 9(k > Ji) (4.7) 

Below we will want to generalize this formula a bit. Towards this end it is worth looking 
back and realizing that the negative sign above arose because we removed Ji boxes from the 
first row - giving a removed tableaux with height 1. Bear this in mind when considering the 
subsequent character formulas we obtain, since we will again obtain 9 functions with a sign 
determined by how many border strip tableau were removed from the first row. 

Using this character formula we have 

(Tr(Z 2j 1)Tr(z2 j 2)Tr(f j 3)) = 2 s xr/Mxr/^) n C * 

R/ihJs iSodd boxes in R 

2 This notation for the Young diagram is listing row lengths, i.e. [k, l j3 ~ k ] has k boxes in its first row and 
then it has J3 — k rows which each have a single box. 
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£^ v ;; (iV-2J 3 + 2£;-2)! 

f (iV + 2fe-2)! 

I (jV_2J 3 + 2fc-2)! 1 ' 



v fc=l k=J 1 + l 

We now consider correlators of four traces. Again, use (I2.15P to obtain 

R/4I-J4 iGodd boxes in i? 

We have J 4 = J\ + J 2 + J 3 and \i is a J 4 cycle while i/ is the product of a Ji cycle, a J 2 cycle 
and a J3 cycle. We know that Xr/4(/-0 wu ^ only be non-zero when i?/4 is a hook. Thus, 
we will consider only Young diagrams R/4 of the form [k, l Ji ~ k ]. Arguing exactly as we did 
above, a simple application of the Murnaghan-Nakayama rule gives 

X[M , <-»]0 1 )X[M ; i- l ]( I/ ) = 6 ( k - Jl ) ~ 6 ( J2 < k < Ji + J 2 ) 

-6(J 3 <k<J 1 + J 3 ) + 9{J 2 + J 3 <k<J A ) (4.10) 

The first term on the right hand side comes from removing both cycles (J 2 ) and (J3) from 
the column. The second term on the right hand side comes from removing cycle (J 2 ) from 
row 1 and cycle (J3) from the column. The third term on the right hand side comes from 
removing cycle (J3) from row 1 and cycle (J 2 ) from the column. The fourth term on the 
right hand side comes from removing both cycles ( J 2 ) and ( J3) from the first row. Thus 

(Tr(Z 2Jl )Tr(Z 2j2 )Tr(Z 2 ^)Tr(Z 2j4 )) = 2 4 £ X R/Mxr/^) J] * 

R/ihJi iGodd boxes in il 

Ji 

= 2 4 J2(°( k ^ J i) - e ( J 2 < A; < Ji + J 2 ) - 9{ J 3 < k < Ji + J 3 ) 



(iV-2J 4 + 2A;-2)! 

* ^ ^3 A \ (jV + 2A; _ 2)! 

16 [Z.- E - E + E J (iV-2J 4 + 2fc-2)! (4U ' 

fe=l fc=J 2 + l fc=J 3 +l fc=J 2 + J 3 + l/ V 7 



Now that we have understood how the signs in the character come about upon applying 
the Murnaghan-Nakayama rule, with a little thought we can obtain the general correlator 
(J n = Ji + J 2 + ■ ■ ■ + J n _i) 

n— 1 / Ji J1+J2 Ji+Jn-i J1+J2+J3 Jn 

(nTr(^)Tr(Z 2 ^)) = T £ - E E + E +-+ E 

4=1 y fc=i fc=j 2 +i fc=j„_i+i fc=j 2 +J 3 +i fc=j 2 +---+j n _i 

(iV + 2A;-2)! 

X (iV-2J n + 2fc-2)! ( 2) 
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This provides a complete generalization of the U(N) result which was derived in jH]. See 

aisoissi sni szi- 

There is a particularly interesting double scaling limit of M = 4 super Yang-Mills 
theory [341 I3"5] defined by 



J 2 

N — > oo and J — >■ oo with — fixed, gyu fixed 



(4.13) 



where J is the number of fields in the trace. In this limit some non-planar diagrams survive, 
leading to a new renormalized genus counting parameter j^. This limit is AdS/CFT dual 
to a pp-wave limit of AdS$ x S 5 , in which the superstring theory can be quantized. Since 
we have computed two and three point correlators we can explore this limit in the SO(N) 
gauge theory. 

In order to extract the double scaling limit of the two and three point correlators we will 
need the following identity 



T{N + Pl + 1) 
T{N-p ) 



(N+ Pl ) 



po+pi 



(N-p -l) 

N Po+Pi+l exp 



11 (N + Pl -l) 



PO+Pl 



E ln ( 1 + 



1=0 



Pi-L 

N ' 



Expanding for large N and summing over I yields 



F(N+ Pl + 1) 



N P0+Pl+1 exp 



r(jv-po) 

Applying it to our SO(N) correlators gives 



^(Pi-Po)(po+Pi + 1) + O(l/N 2 ) 



, m lB » lsm / ^27^^ T(N + 2(J — k) + 1) 

Tr (Z 2J )Tr(Z 2J )) = 4 > K — — J — 



4iV 2J exp I 
k=i ^ 



( J{2J -Ak + 1) 



N 



4JiV^ Sinh - 



2.P 

N 



and 



(Tr (Z 2Js )Tr (Z 2Jl )Tr (Z 2J *)) = 8 



Ji 



E- E 

v fc=i fc=j 2 +i, 



r(iy + 2(j 3 -fc) + i; 
r(Ar-2fc + i) 



9 , Ja sinh %^ sinh 2%Za s i n h 2^/3 sinh 2^ 
8 x 4iV 2j3 e^ ^ ^ 4 x 2J,iV 2j3 * 



e n — 1 



A/ 



Compare these results with U(N) correlators in the double scaling limit 
(Tr (Z 2J )Tr (Z 2J )) = £ ^E+jL „ 2 JiV 2J ^ ^ 



-? T(N-2J + k) 
13 



2J£ 

AT 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



<TY(Z"3)TY(Z" W Z 2A )>=( E -E) rw-^*l» 

\fc=2J 2 +l k=lj ^ 3 ' 

~ 2A^e^ (coth ^ - 1) sinh ^ sinh 2 J2 J3 



A" A" iV 

9 . sinh %^ sinh 
~ 2J 3 A^ 2J3 L ^- (4.19) 

£3 

A' 

It is clear that non-planar unoriented diagrams in SO(N) gauge theory do not survive this 
limit. 

Our result is similar to earlier results obtained for the double scaling limit of the matrix 
model relevant for the c = 1 string [49~]. From this point of view, the ribbon graphs of the 
matrix model are identified as a triangulation of the string worldsheet. These double scaling 
limits take iV — > 00 simultaneously with the world sheet continuum limit in such a way that 
the string coupling is held finite, so that sums over continuum surfaces of any topology are 
captured. This limit for antisymmetric matrices has been discussed in [50]. In this double 
scaling regime too, only orientable surfaces survive. 



5 Counting 

When constructing gauge invariant operators, it is possible to contract indices using any in- 
variant tensors. For SO(N) there are two such tensors: the Kronecker delta (which contracts 
pairs of indices) and the e %ll2 "" lN tensor. In [TJ, taking both invariant tensors into account, 
we made a precise conjecture for the basis that can be constructed. In this section we would 
like to count the number of operators we proposed and thereby verify that it is indeed a 
complete set. We will focus on the case that iV is even. We will end this section with a few 
comments on the odd iV case. 

To start, recall the discussion of [TJ: operators of the form ( II .ip are the complete set of 
operators that can be built without using e ilt2 "' iN . They correspond to the set of operators 
that can be written as a product of traces of even powers of Z. Further, they are labeled by 
Young diagrams with an even number of boxes in each column and row. As explained in [1] , 
these Young diagrams R built using 2n boxes with the restriction that R has not more than 
A" rows, can be indexed by partitions of | that have no more than y parts. Consequently, 
we can write the partition function for the number of operators of the form ( II. ip as 

8=1 

The coefficient of the x p in the expansion of Fi(x) tells us how many operators can be 
constructed using p Z fields. Next we need to count the number of operators (11.21) built using 
e %\ii-% N _ These operators all have a dimension > y. They are constructed by contracting 
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the indices of y Z fields with e* 1 * 2 "' %N and contracting the remaining indices in pairs. These 
can also be constructed in the form (jl.lj) . with labels i? that consist of a single column of N 
boxes with a second Young diagram stacked to the right. To get a non-zero operator, this 
second Young diagram must again have an even number of columns and rows. Clearly then, 
the partition function for the operators constructed using one e llt2 '" %N is 



JV 
2 

JV 



F 2 (x) = x^F 1 (x) =^\[-^Y i (5-2) 

i=i 

Operators constructed using any even number of e Hl2 "" tN s lead to Young diagrams that have 
both an even number of rows and columns, and hence are included in (15. ip . Similarly, 
operators constructed using any odd number of e lin "' lN s are included in (15. 2p . Consequently, 
the complete partition function for the operators in our basis is 



JV 
2 

„ , N , JV 

Fix) = (1 + x~ 



II. < 5 - 3 ' 

i=l 



The partition for free Yang Mills theory on a compact space has been computed in 
The result is 

/do r E 

[dO] x nE Xsy^ (R) (0) = / [dO]e^ ^^(° m ) (5.4) 

n=0 J 

Here we take R to be the adjoint representation (since our field Z transforms in the adjoint) 
and Sym n (i?) is the representation obtained by taking the symmetric product of n copies 
of the adjoint. Set E = 1 and then expand to get a polynomial in x. The coefficient of x n 
counts the number of operators that can be built using n Zs. Consequently, if our basis is 
complete, we should find G(x) = F(x). To evaluate (15. 4p we need the adjoint character (we 
focus on SO(2n)) 

Xr(x) = (XjXj + x^Xj + XjXj 1 + x^xj 1 ) + n (5.5) 

l<i<j<n 



and the integration measure valid for any symmetric /(x), x = (xi, .... x^,)[5T] 

/, ' do|/(x) = f nJ|: A(x+x ~ 1)2/(x) (5 ' 6) 

^SO(2n) Z n - JT n =1 ^TT^Xj 

Here T n = S 1 x S 1 x • • ■ x S 1 is the unit torus. Using this we find 
G(x) = [ [dO] e £~=ii^xB(o») 



'SO(2n) 

1 f W dx 



2 n ~ 1 n\ At 11 27rix 



1 1 9lT7V • V 



-1\2 
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E™ = i^(Ei< l < J <„(^r^+^ m ^+^r^ ra +^" m ^ m )+ n ) 

-1\2 



xe 

1 f TTJ^A(x + x- 1 ) i 
J Tn f = \2mx 3 1 J 

nil 1 1 1 , , 
— r-r rr — r (5.7) 

l<i<j<n J L J J 1 J 

Some straight forward algebra shows that 



2«- 1 n! 



>/J " j=l J l<i<i<n v y 



where 



*^ 1 . . 1 / >-» / >-» 1 rr> 1 / >-» net / >-» 1 1 rp rp 1 / >-» 1 

We have not managed to compute these integrals for general n. We have however verified 
that (the integrals over Xi below all run over the unit circle; we have taken x < 1) 

dxt f dx 2 1 

F12 = 73— — (5.10) 



4 J 2nix\ J 2nix 2 (1 + x) 

so that for iV = 2n = 4 we have 



G(x) = (l+x 2 ) r ^ I ^ I (5.11) 



and 

1 /" d Xl r dx 2 r dx 3 = 1 1 1 

24 J 2m Xl J 2nix 2 J 2mx 3 12 13 23 (1 + x f \ + x 2 x + x + ^2 ^ • ^ 

so that for iV = 2n = 6 we have 

G{x) = (1 + x 3 ) . 1 21 1 41 1 6 (5-13) 
1 — x z 1 — x 4 1 — X D 

This proves that for iV = 4, 6, G(x) = F(x) which provides strong support that we have 
indeed constructed a basis. Given the results of this section, we conjecture that the partition 
function for SO(N), with N even is given by (15. 3p . Our results also lead us to conjecture 
that for SO(N), with N odd, the partition function is 

N-l 

= fl r=b^ (5 - 14) 

i=i 
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6 Link to Free Fermions 



The dynamics of gauge invariant operators of a single Hermitian matrix model reduces to 
eigenvalue dynamics[37]. If one restricts to the gauge invariant and purely holomorphic or 
antiholomorphic observables of complex matrix models, one can again reduce to eigenvalue 
dynamics |52j. This eigenvalue dynamics can then be mapped to the dynamics of iV free 
fermions in an external potential[2]. For the free complex matrix model, the basis provided 
by the Schur polynomials has a particularly close relationship to free fermion dynamics: the 
Schur polynomials can be mapped to free fermion wave functions j2]. In this section our goal 
is to argue that our operators (\3.7\i also describe free fermion wave functions. 

Before considering the eigenvalue dynamics it is useful to recall a few facts about the 
wave functions of the single particle Hamiltonian 

d d 

which describes a single particle moving in a harmonic oscillator potential, in two dimensions. 
The ground state wave function (not normalized) is 

M*,z) = e~ 22 " (6.2) 

Further, at any given energy level, the state with largest angular momentum (again, not 
normalized) is given by 

ipi(z,z) = z l e' z ~ z (6.3) 

The parity of this wave function is given by (—1)' so that (for example) if I is even we have 
an even parity state. These wave functions will play a prominent role below. 

To start, return to the two point function (13.81) and rewrite it in terms of eigenvalues 



<Wtf,odd = (xUZ)Xs(Z)) 



N 



[dzdz}\A(z)\\ R (z)x s (z)e- 2 ^ z ^ (6.4) 

where the required Jacobian A(z) = Yli<i<j<Z-( z i — z ]) * s discussed in Appendix O Conse- 
quently, if we introduce the wave functions 

N 

^ B {z,z) = A(z) X ° R (z)e-^ z ^ (6.5) 

the two point equation (16 .4p becomes the statement that these wave functions are orthogonal. 
To find the interpretation of these wave functions, we need to unpack some of the details 
of the expression ()3.7|) . Study operators built from n fields. Recall that n must be even. 
We will focus on SO(N) gauge theory with N even. Denote the eigenvalues of Z by ±iz{. 
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Introduce the matrix 



M 



z% 












z 2 



(6.6) 



Using the result (IC.lOp from Appendix [UJ we have 

(-l)f X5/4 (M) 



X°s(z) 



(6.7) 



where xs/a{M) is nothing but the Schur polynomial. The Young diagram S can also be 
taken as a partition of ~ with parts s\, S2, sn equal to the row lengths of S. Recall that 

z 2 

the Schur polynomial obeys the identity 



xs{yi,y2,-~,y*] 



I j (Vi - Vj) x det 



l<i<i<^ 



V\ V2 



Vi 



V2 



si- 

2/jv 

2 

S2- 
2 



JV 
' _> " 



!Jl 



Using this identity we find 



-l)2det 



ri+iV-2 ri+N-2 
z l z 2 

r 2 +N-4 ^r 2 +N-4 



Z 



Z 



I- N 
„ "3" 



v T 



n+AT-2 

^ JV 

r 2 +AT-4 

* JV 



''/V 

^ JV 

_> 



2/jv 

2 



(6.8) 



' Xj-i — l -^i-^i 



(6.9) 



The interpretation of this wave function is clear. It is the Slater determinant of ^ single 
particle wave functions, and consequently is the wave function of y fermions. These fermions 
are not interacting (apart from the usual Fermi statistics) but are moving in an external 
harmonic oscillator potential. Each particle has maximum angular momentum for its energy 
and is in a parity even state (recall that the number of boxes in each row of R is even). 

Next consider the wave functions ipr(z, z) that corresponds to the operators Qr(Z) given 
in (11.21) . Since the approach to this problem is identical to what we have done above, 
it should be no surprise that the wave functions <Pr{z, z) take exactly the same form as 
( 16. 9p . The key difference is that now every row has at least one box and all rows have an odd 
number of boxes. Thus, we again have y fermions, each particle again has maximum angular 
momentum for its energy, but now each is in a parity odd state. This result should not be 
surprising. We can describe these fermion wave functions in terms of a density in phase 
space, indicating which states are filled. For some nice examples, see [33]. The operators 
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Qr(Z) include a factor of the Pfaffian. In terms of the fermion phase space density, the 
Pfaffian creates a small hole at the origin of phase space, pushing all fermions up by one 
level so that each particle lands up moving from a parity even state into a parity odd state. 

There is a simple generalization of this result to SO(N) with iV odd. The two point 
function (13. 8 p becomes 

<W*,odd = (x°r(Z)xs(Z)) 

/JV-1 
[dzdz]\A(z)\\ R (z) X ° s (z)e~ 2 ^ '** (6.10) 

JV-1 

where now A(z) = Ylk=i z k Yli<i<j<^-( z i ~ z j)- This Jacobian is discussed in Appendix O 
A very similar discussion shows that the corresponding wave functions are unchanged. This 
follows because in the formula (16. 9 p we need to replace A" — > N — 1 and we need to multiply 

JV-1 

by the extra factor Yik=i z k- Multiplication by the extra factor can be accomplished by 
sending rj — > Tj + 1. The combined transformation A" — > N — 1 and r j — > r j + 1 leaves (16. 9p 
unchanged. Since A" is odd, the fermions are now all in a parity odd state. 

The identification of our operators with free fermion wave functions has an immediate 
application. There is a class of \ -BPS bubbling geometries which arise as orientifolds of 
type IIB string theory. In [53] these geometries have been put into correspondence with 
free fermions (see also [54]). Using this dictionary, we are able to give the operators dual 
to certain backgrounds with these geometries. The description of [53] is in terms of free 
fermions moving in a harmonic oscillator potential on the half-line. The description we have 
constructed above is in terms of the holomorphic sector of two-dimensional free fermions (i.e. 
2+1 dimensional) moving in a spherically symmetric harmonic oscillator potential. The two 
are related by truncating the system of two dimensional fermions (analogous to projecting 
to the lowest Landau level in the quantum Hall effect). See [2J [55] for a clear discussion of 
the relevant truncation in the context of the U(N) gauge theory. 



7 Sp(N) gauge theory 

There is a close relation between invariants of representations of SO(N) and representations 
of Sp(N) obtained by exchanging symmetrizations and antisymmetrizations and replacing 
A" by — N. In this section we will argue that the finite A" physics of SO(N) and Sp(N) 
gauge theory are related in exactly the same way. 

In Sp(N) gauge theory, we have the following identity obeyed by fields in the adjoint (X 
is an N x A" matrix with N even) 

x T = -nxn T = nxn (7.i) 
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where the N x N matrix Q is given by (1^ is the d dimensional identity matrix) 

f In 

n=\ ^ 
\-H 

The free field two point function is 

(Xij X k {) = 5u5jk — Qik^ji (7.3) 

Note that Q clearly satisfies the following relations 

tt 2 = -1 N , = tt T = -a (7.4) 

A clear consequence is that only traces with an even number of fields survive 

Tr (X J ) = Tr ((X J ) T ) = (fi 2 ) J Tr (X J ) = (-1) J Tr (X J ) (7.5) 

Now, consider a two matrix model with real fields X and Y . Combine these real fields into 
the complex combinations 

Z ab = X ab + iY ab^ ^ = x ba + iy ba (j 

The free field Sp(N) two point functions are given by (Z^Z H ) = 0, (ZijZki) = and 

(z^z kl ) = sisj + (ny l (n- 1 y k (7.7) 

This looks very similar to the two point functions of the SO(N) gauge theory (see figure H]). 
Consider the computation of (T°* u ) where 

in the free field theory. Our goal is to relate the computation of this correlator in the 
Sp(N) theory with the computation of the same correlator, but in the SO(N) theory. The 
contribution coming from keeping only the first term (= 5 l k 5l) in each Wick contraction (see 
figure [2]) is identical for the two. This term has the form N p . The contributions coming 
when we start to include the second term (= (fi)J(fi -1 )], for Sp(N) and = — Sftl for SO(N)) 
in the Wick contractions can differ in sign. Contributions coming from including any odd 
number of these second terms will differ in sign between the Sp(N) and SO(N) theories, 
and will have the form cN p ~ 2l ~ l for some integer i and c a constant determined by the 
number of ribbon graphs summed. Contributions coming from including any even number 
of these second terms will have the same sign in the Sp(N) and SO(N) theories, and will 
have the form cN p ~ 2k for some integer k. Thus, when p is even the computation of (17. 8p in 
the Sp(N) theory gives the same answer as the answer for the SO(N) theory after we take 
iV — > —N. When p is odd the answer for the Sp(N) theory follows by taking iV — > —N 
in the answer for the SO(N) theory and multiplying by —1. The power p of the leading 
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term is determined rather simply by the sign of the permutations /3 and r. To compute the 
sign of a permutation decompose it into a product of transpositions. This decomposition 
is not unique. The sign of the permutation sgn(a) = (— l) m where m is the number of 
transpositions in the product, sgn(a) is well defined, i.e. it does not depend on the specific 
decomposition of a into transpositions. Our final result is 



(T^SPiN) = sgn(/3)sgn(r) 



N^-N 



This result will play an important role below. 



(7.9) 




+ 




(A) 



(B) 



Figure 2: A comparison of the free field two point function of the SO(N) gauge theory 
(shown in (A) above) and of the Sp(N) gauge theory (shown in (B) above). 



Introduce the operators 

o R (z) 1 



(2n)! 



«l*2*2*l"'«n-l*n«nin-l 
jlj2"'hn-lhn 



. . . ^2n-lJ2n 



O r (Z) 



(2n 



_ v i \ hh-h^-ih-n, y y 

\\ / j X.R\ U ) U hi2i2h— in-lininin-l ^hh 



J2n-U2n 



(7.10) 



(7.H) 



0"£S2r, 



To compute the two point function of these operators we will be using (17. 9p . For this reason 
we will need to spell out the gauge group we are using to compute the correlator. We will 
need one more result from group theory. Recall that to get the conjugate (or transpose) R c 
of a Young diagram R we need to swap rows and columns. For example 



R 



R c 



The characters of R and R c are related by 



(7.12) 



(7.13) 
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From the point of view of the projectors used in defining R} taking the conjugate of R 
corresponds to swapping symmetrization and antisymmetrization of indices. We are now 
ready to compute the correlation functions of the R . For N even, R h 2n and S h 2m we 
find 



— — r r 

(2n)\ (2m)! ^ 



(0 R (Z)O s (Z)) Sp(N) = — — — — y y 

I i\ilili\-in-\in^n-\ yjxh yh-n-\h-n T nn---nm-\nm y y \ 

\ hh-hn-ihn i1i2i2i1-im.-1imim.im-1 n hh n hm-ihmi 'Sp(N) 



(2n)\ (2m)! *-f ^ 



XR(cr)xs(r)sgn(a)sgn(T) 



v / »l*2*2tl—tn-l»n»n*n-l yjlh yjiji yj2n-lj2n r JU2---J2m-lJ2m 7 7 \ 

\ jlh-hn-lh-n ' ' ' ili2i2il-im-limimim-l hh ' ' ' Zj j2m-lj2m/ SO(N) 

v 1 v y <xes 2n re5 2m 

; ili 2 j 2 il---j„_li„i„i n _l r^ja . . . yj2n-lj2n T h32---j2m-lj2 
\ jlh—hn—lhn " ' ' ili2i2h—im-l 

= [(0 R o(Z)Osc(Z)) so{N) 



U2m 7 7 \ 

limimim-i hh ' ' ' ^ ' hm-lhm I SO(N) 



(7.14) 



Using the results of [I] we now immediately fincjfl 



(0 R (Z)O s (Z))s p{ N) = S RS 2 n 
It is straight forward to verify that 

n 



l R/4 

d R 



n 



tGodd boxes in R c 



(7.15) 



iGodd boxes in R c 



n 



(7.16) 



iSeven boxes in R 



so that we finally obtain 



(o R (z)o s (z)) Sp(N) = s RS r 



lR/4 

d R 



n 



(7.17) 



igcvcn boxes in R 



Recall that the even boxes will occupy every second row, including the top most row. As an 
example we have filled the even boxes below 



(7.18) 



We can now repeat many of the arguments we developed for the SO(N) theory. Since 
many of the details are basically the same, we will for the most part simply quote the results. 

3 Recall that as symmetric group representations we have da = cIrc . We have used this when we wrote 

(EE). 
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Correlations functions are given by 

(Tr Tr (v{Z 2 )^)) Sp{N) = 2 1 ^ 1 ^ £ Xr/Mxr/^) R Q (7.19) 

R/A\-n/2 iecvcn in R 

Using this result we can again give a simpler form for our operators 

X°s(Z) = 7^7 E 2-< M X^)Trv^M^T f ) (7-20) 

The two point function of these operators is 

(x° R (Z)x s(Z))s P (N) = 5ns ]J c % (7.21) 

iSeven boxes in S 

These operators again enjoy a simple product rule (S/4 h R/4 h ^) 

X° 5 (^)xM = 9R/AS/AT/4J&W (7-22) 

T/4 |_ 211+22 

where gR/is/iT/4 is the Littlewood- Richardson coefficient. 

Using (I7.19P we can again give a formula for general extremal correlation functions 

(Tr (Z 2J ")Tr (Z 2Jl )...Tr (Z 2 ^- 1 )) Sp(N) = 



fc = l fc = J 2 + ...Jn-l + l, 



where J ra = Ji- m particular, for two point functions we find 

We have also considered the double scaling limit of the two-point functions in the Sp(N) 
gauge theory. We find 



(Tr (Z 2J ) Tr (Z 2J )) ~ 4iV 2J ^exp 



fc=i 



^(4*-2J-l) 



J(2J+3) ,7(2.7-3) . . n r2 

, e jv — e ^ . (-j r smh 



4 jV 2J - _^ 4JiV 2J 9I2 N (7.25) 

— 1 



2J2 
N 



Similarly for three-point functions 

r(JV + 2fc) 



/ Ji Js \ 

(Tr (Z 2j3 )Tr (Z 2Jl )Tr (Z 2Ja )) =-8 J- ^ 

\fc=i k=j 2 +i/ 



T(N-2(J 3 -k)) 



32iV 2j3 ^ . 2J 1 J 3 . 2J 2 J 3 _ , AT2J3 sinh ^ sinh 

-e * smh _ smh _ ~ 8 J 3 iV ^ (7.26) 

23 



Since this is the same as the SO(N) result, we again see that in this double scaling limit 
only orientable higher genus surfaces contribute. 

Our operators can again be related to free fermion wave functions. Indeed, start from 



rsjr ,cvcn 



(x°r(Z)x°s(Z)) 



[dzdz}\A(z)\\° R (z)xUz)e 



2 — l 



(7.27) 



where A(z) = Uk=i z k IIi 



zf — z 2 -). Consequently, 



V 



if> R (z,z) = A{z)x R {z)e-^ z ^ 
Arguing as we did for the SO(N) theory, we find 



i[>r(z,z) = det 



z 



n+N-l ri+JV-l 



Z 



Z 



r 2 +N-3 „r 2 +iV-3 



rjv+1 
v "2" 



Zn 



r N +l 
v "2" 



n+N-l 

* JV 



r jv +1 

Z 

2 



(7.28) 



(7.29) 



This is the Slater determinant of y single particle wave functions i.e. the wave function 
of y fermions. These fermions are not interacting but are moving in an external harmonic 
oscillator potential. Each particle has maximum angular momentum for its energy and is in 
a parity odd state (recall that the number of boxes in each row of R is even and iV is even). 
This completes our discussion of the Sp(N) gauge theory. 



8 Discussion 

There are a number of interesting questions that could be pursued at this point. In this 
section we will list some of them. 

Using the results of this article, it is possible to study the three-point function of two giant 
gravitons and one point like graviton in the SO(N) gauge theory. These correlators can also 
be computed using methods of semiclassical string theory, using a Born-Infeld description of 
the giant graviton. Computations of this type[56J \57\ EH1 EH] show a perfect match between 
correlators computed in the U(N) gauge theory and three point functions computed using 
semiclassical string theory. Computations along these lines will further test the gauge theory 
results we have obtained. 

It has been possible to construct a very detailed holographic map between ^-BPS bub- 
bling geometries [54J and M = 4 super Yang-Mills theory with U(N) gauge grouppO]. The 
connection between the |-BPS sector of Af = 4 super Yang-Mills theory and the dynamics 
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of free fermionsj^l EI] was a key insight that determined much of the holographic dictionary. 
Given the close connection to free fermions developed above, it is natural to develop the 
holography of bubbling orientifolds of type IIB string theory. 

Finally, it would be natural to study the spectrum of anomalous dimensions in non-planar 
large N limits of the SO(N) and Sp(N) gauge theories. New integrable sectors of the U(N) 
gauge theory have been found in these limits [T4" | fT5l fT6] . Is integrability also present in large 
TV but non-planar limits of SO(N) and Sp(N) gauge theory? 

We hope to return to these issues. 
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A Simplifying the SO(N) basis 

In this Appendix we will give an alternative derivation of ( 13. 6p . Our starting point is 

^4E x R w)c;x.., 2n witt;^ 11 ' 2 • ■ ■ zi2n - li2n (a- 1 ) 



2n! 

Recall our shorthand notation for indices 



c^ = c^ , , (f3) T = (l3) nn '" J2n , Z 1 = Z hi2 ■ ■ ■ Z i2n - li2n . (A.2) 



First, the contractor Cj 4m can be written as 



5- 6 • 



8 8 8 
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The second line follows because cx4 M (r) = r + 1 for r odd. Consequently, (lA.ip becomes 

ot = hY xr^ctKz 1 



2n\ 



2n! 



«-y E X/?(/3)4 1 fc 2 ---4 2n _ 1 fc 2n (/3^4 M )f^ / 



2n! 

/3e5 2n 



-y XniPvii^kika ■■•4 2n _ 1 fc 2 „(/3)f Z 7 - 



2n! 

/?e5 2 



(A.4) 



Since Z is antisymmetric, for any 77 G S^fS^] we have 

^(/) = ^ ( i)S(2) . . . ^ (2 „-i)S( 2 „) = sga^z 1 , (A.5) 

It is now straight forward to see that 

4 1 fe 2 ---^ 2n _ 1 fc 2 „(^)f^ / = S klha ---6 klln _ lhan Z r, ~ lfi ~ l6 ~ 1 ^ 

= S klk2 ■ ■ ■ 5 k2n _ lk2n (Pri)j Z T sga(C) 

= S klk2 ---6 k2n _ lk2n (/3)fz I sga(0. (A.6) 

Consequently 

Or" = ^ XR(^4^k 1 k 2 ---S k2n ^ k2n (f3)fz I . 

' /3eS 2n 

' K '' /3es 2 „ ^ v eS n [S 2 ] 



1 1 



2n\ (2 n n\) 2 

Now, 0A.6P implies that the expression within brackets defines a function on the double 
coset S'rJSy \ S2nl S n [S^[, i.e. it takes the same value for different a G S^n that represent the 
same double coset element. Thus, we can trade the sum over (3 G S211 for a sum over coset 
representatives which is a sum over partitions 

0°r = 2~j Y Y ~ X«(^4^o-^ 1 )sgn(04ifc2---^2„-ife2„(^)f^ / - 

= ^ Y — 2KM)+ ' H Xi?/4(M)Xi ? /4(^)^ lfc2 ---4 2B _ lfan (/3 4 ,)^ 7 , (A.7) 
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To obtain the last line above we have used the mathematical identity 

^ E X«(^<)sgn(0 = 2 ,W+iW X«/4W»/4M. (A.8) 

Thus, our operators become 

U R - Z 1 >. 1 XR/A V ) d klk2 ■ ■ ■°k2n-lk2n{P±l<) I Z 

7 — /„ Z 4v 



vhn/2 



= 2 ' Mr ^7^yy E Z^'W^W^T*)- (A.9) 

For the special case /i = (l n//2 ) we find 

Oa = 2^ d -P T ^- E 2^> W (a)Tt(a(*)°*). (A.10) 

which completes the demonstration. 

B Correlation Functions 

To test the result ( 12.15)) we have studied a number of correlation functions using Mathematica 
and analytic techniques. Consider first correlators of the form (Tr(Z 2 ) p Tr(Z 2 ) p ). Introduce 
the notation 

A p = (Tr(Z 2 ) p Tr(Z 2 ) p ) (B.l) 

By studying Wick contractions (or equivalently Schwinger- Dyson equations) it is not hard 
to obtain the following recursion relation 

A p = (16p(p - 1) + ApN(N - 1))A P _! (B.2) 

Since Aq = 1 we can easily generate explicit expression for A p . For example 

A 5 = 122880iV 10 - 614400iV 9 + 6144000iV 8 - 20889600iV 7 + 98918400iV 6 - 226222080iV 5 
+60456960(W 4 - 855244800iV 3 + 1148190720iV 2 - 754974720iV (B.3) 

Introduce the short hand 

II c * = fn d ( B - 4 ) 

iSodd boxes in R 
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Our formula ( I2.15P says 



A, 



,10 



fl 



odd 



odd 



+ 25/ r 



odd 



+ 36/ r 



odd 



+ 25/ r 



odd 



+ 16/r 



odd 



V 

which is indeed correct. 

Next, consider correlators of the form (Tr(X 2 ) p Tr(X 4 )Tr(X 2 ) p+2 ). Introduce the notation 

B p = (Tr(X 2 ) p Tr(X 4 )Tr(X 2 ) p+2 ) (B.6) 

Again by studying Wick contractions (or equivalently Schwinger-Dyson equations) it is not 
hard to obtain the following recursion relation 

B p = 16(2p + 4)(2p + 2)(2p)(2p - 2)5 P _ 2 + 2(2p + 4)(JV - 
+ [4N 2 {N - l)(2p + 4)(2p + 2) + (2p + 4)(2p + 2)(2p)8(3X - l)]A p (B.7) 

We now easily find, for example, 

B 3 = 245760iV 9 - 1105920X 8 + 10813440iV 7 - 32686080X 6 + 142786560X 5 

-277831680X 4 + 586383360X 3 - 617349120iV 2 + 188743680X (B.8) 

Our formula (12.151) says 

/ \ 



/r 



odd 



+ 



podd 



+ 5/r 



odd 



+ - 5/r 



odd 



podd 



fl 



odd 



B.9) 



V 

which is indeed correct. 

Finally, we have also performed a complete check in Mathematica for the full set of 
operators that can be built using 8 fields. The Mathematica and analytic results from (I2.15P 
are again in complete agreement. A study of the Sp(N) coorelators has also been performed 
to confirm (I7.19P using Mathematica. 
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C Jacobians 

For single matrix models, since we are interested in the dynamics of gauge invariant observ- 
ables, we can employ an eigenvalue description. The idea is to write the path integral as 
an integral over the eigenvalues and some angles. We then integrate out the angles. The 
resulting measure is nontrivial. With a slight abuse of language, we refer to the measure as 
a Jacobian. In this Appendix we would like to compute the Jacobian for a single complex 
SO(N) matrix and for a single complex Sp(N) matrix. To the best of our knowledge, these 
are new results. 

The approach we employ for determining the measure, is first to compute it for a real 
matrix and then use this result to guess the answer for the complex matrix. We check this 
guess by verifying that we get the correct answer for any correlation function we compute. 
To determine the Jacobian for a real (or Hermitian) matrix, we require that the Schwinger- 
Dyson equations in the original (matrix) variables agree with the Schwinger- Dyson equations 
in the eigenvalue variables. This implies a differential equation for the Jacobian which we 
solve. Although the Jacobians we compute in this way are all known, we have given our 
derivation since it seems to be new and is simpler than existing derivations. To start, we 
illustrate the method with the U(N) matrix model and then move on to SO(N) and Sp(N). 

The use of Schwinger- Dyson equations to determine a Jacobian in this way was pioneered 
in collective field theory [62]. See (63] for applications to multi matrix models and [M] for 
applications to vector models. 



The main goal of this subsection is to illustrate how we compute the Jacobian using Schwinger- 
Dyson equations and also to illustrate the close connection between the Jacobian for the real 
(or Hermitian) matrix and the complex version. 

Consider the matrix model for a single matrix X living in the Lie algebra u(iV), i.e. X 
is Hermitian. Rewrite the Schwinger-Dyson equation 



C.l U(N) Matrix Models 








(C.l) 



in terms of eigenvalue variables to obtain 
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where J(A) is the Jacobian we want to determine. After performing the sum over r we have 



N N N 



0= f[dx]j(x) 2 T^- + (--i)E A r 2 -E A n e ~ IE£lA? ( c - 3 ) 



i,j=\i^j J i=l i=l 

Now, work directly in the eigenvalue variables 

N 



= fmY,w.(*r 1 JWe~ i ^) 
= W ( A)(f xr^§^ + (»-i) eat" -Eat) e - E - A? to*) 

Comparing (1C.3I) and the second line of flC.4j) we learn that 



dlogJ(A) A 



_ 2 E T~3T~ ( c - 5 ) 

This clearly implies that 



OXi ' Aj — A^ 



which is indeed the correct result. 

The argument given above is not quite rigorous. Indeed, to really prove that we will 
get the correct answer for any gauge invariant observable in the theory, we should have 
considered the Schwinger-Dyson equation which follows from 

J IJ 

(C.7) 

with S = |Tr(X 2 ) and with J 7 any gauge invariant observable in the theory. This gives the 
complete set of Schwinger-Dyson equations of the theory. This much more general argument 
continues to imply ( 1C.5|) . Further, we could also have included interactions (i.e. considered 
"actions" S with higher than quadratic terms); one again finds (IC.5j) . In what follows, for 
the sake of simplicity, we will continue to work in the free theory and simply set J 7 to one. 
The reader should bear in mind that adding interactions and/or considering the most general 
T do not change our results. 

We will need a generalization of the above result. Consider a matrix model of two 
Hermitian matrices X and Y. Collect these into the complex matrix Z = (X + iY) and 
its Hermitian conjugate Z*. We are interested in correlation functions of traces of Zs or 
traces of Z^s, but not of traces in which both Z and Z^ appear in the same trace. Denote 
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the eigenvalues of Z by z\. The eigenvalues of are then 2^. The Jacobian in this case 
is I 



J = A(z)A(z) 



The relation between (1C.6P and (IC.8P is striking; we will find a similar relation between the 
Jacobians for the real and complex SO(N) and Sp(N) theories below. 

C.2 SO(N) Matrix Models 

We will now consider the matrix model relevant for SO(N) gauge theory. We consider the 
case that X is even and start with a single real N x N antisymmetric matrix X. In this 
case, we can, with a unitary transformation, bring X into a block diagonal form 



X 



xx 

-xi 











x 2 

-x 2 



(C.9) 



Using this explicit form for the matrix X we easily find Tr(X 2J+1 ) = and 

N_ 

2 

Tr(X 2J ) = 2(-l) 3 x2 i J ( C - 10 ) 



»=i 



Since X is an antisymmetric matrix, elements above and below the diagonal are related by 
a sign. Consequently 



dX. 



dX, 



5ik5ji — 5u5jk 



(C.ll) 



To get some practice with this derivative and since we need it in what follows, consider 

2J-2 



d 



dX 4 



<x 2J -% 



dXj a 

2J-2 



dX ki v 2J- r -2- 



'^^(X r ) ik (5ik8ji — 5u5kj)(X 2J ' 2 )ij 

r=0 
2J-2 

[Tr(X r )Tr(X 2J - r - 2 ) - Tr((X T ) r X 2J ~ r - 2 )] 

r=0 

J-l 2J-2 

^Tr(X 2r )Tr(X 2J - 2r - 2 ) - ^ (-l) r Tr(X 2J " 2 ) 



r=0 



r=0 
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J-l 

^Tr(X 2r )Tr(X 2J - 2r - 2 ) - Tr(X 2J " 2 ) 

r=0 



Now, consider the Schwinger- Dyson equation 



''-%ei^ x "' 



which implies 

j-i 



(^Tr(X 2r )Tr(X 2J " 2 ^ 2 ) - Tr(X 2J ~ 2 )) = -2(Tr(X 2J )) 



r=0 



Writing this in terms of eigenvalues gives 

j ^ N_ JV JV 

(4(-l) J+1 E E 'f •'7 / " 2 - 2(-l) J+1 X:^- 2 ) = -4(-l)'£> 

r=0 i,j=l 1=1 i=l 

or, after summing over r 

AT JV JV 

2 2J 2 2 

< 8 E ^+ 2 ( 2 ^- 1 )E^ 2 > = 4 <E* 2J > 

i,j=l,i^j 1 3 i=l i=l 

Now, working in terms of the eigenvalue variables we have 

JV 



which becomes 

/JV 2 / 

i=l ^ 

Comparing (1C.37|) and (1(1390 we have 



xf- l ^- + {2J-l)x 2 l J - 2 -2xf 



d log J 4xi 



which is solved by 



2 



m= n (^-^ 2 ) 2 = (A(x 2 )) 2 

i<j=l 
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This is the correct answer. 

We now consider a model of two real SO(N) matrices X and Y. Introduce the complex 
matrix Z = (X + iY) and its Hermitian conjugate Z* . Denote the eigenvalues of Z by 
Zj. It follows that the eigenvalues of Z^ are given by complex conjugation as Zi. We are 
interested in computing correlation functions of multi trace operators that have only Zs or 
Z^s in each trace. Based on the relation between (IC.6p and ( 10.8j) . we conjecture that the 
Jacobian is given by 



J(z,z) = A(z 2 )A(z 2 ) 



(C.21) 



To prove that this guess is correct, we will compute two point functions of our operators 
Xr(Z) and show that they agree perfectly with our result ( 13. 8 j) computed directly in the ma- 
trix model. Since the most general multi trace operators that have only Zs or Z^s appearing 
in each trace can be written as a linear combination of the X°r{Z) and XhC^ )? this proves 
that (IC.2ip is the correct Jacobian. 

Correlation functions in the eigenvalue basis are computed using 



(•••)= N \ [dzdz]J(z,z)e 



(C.22) 



We will fix the normalization by requiring that (1) = 1. To perform this computation, note 
that we can write 



a(z 2 ) = n^ 2 -*?) 
i<j=i 



det 



y N-2 
'1 

y N-4 
'1 



4 



N-2 
z 2 

~iV-4 



r 2 
-2 



■ni2-iN N-2h „N-2i 2 



T z 



-9 



N-2 

" N_ 
2 

JV— 4 

6 jy_ 
2 



z 2 

<6 JV 



N-2i N 
Z jv 

2 



(C.23) 



We will make frequent use of the identity 



dzdz z n z m e 



7iS mn n\ 



(C.24) 



It is now straight forward to verify that 



Af 



(C.25) 



33 



Next, using the identity 



A(z 2 ) X ° R (Z) = ^"-n . .. Z y +N ~ 2 " (C.26) 



2 



we easily find 

JV 

(XUZ)X° S (Z')) = U ^ + N ~ 2l)l = 6 RS f RMd (C.27) 

rLLi(iv-2z)! 

in perfect agreement with ( 13.81) . A simple generalization of this argument shows that we 
correctly reproduce the correlators for odd N too. For the odd N case the Jacobian is 

J{z,z) = A(z)A(z) (C.28) 

where 

JV-l 

fc=i i<i<j<^i 

C.3 Sp(N) Matrix Models 

Finally, consider the matrix model relevant for Sp(N) gauge theory. In this case, N is even. 
Start with a single real N x N matrix X in the Lie algebra of Sp(N) which hence obeys 

X T = ttXtt (C.30) 

It is easy to see that 

= det(Al - X) = det(Al - X T ) = det(Al - ttXtt) = det(Al + X) (C.31) 

so that the eigenvalues of X come in pairs ±Xj, i = 1,2, Consequently, we find 
Tr(X 2J+1 ) = and 

N_ 

2 

Tr(X 2J ) = 2^a; 2J (C.32) 

8=1 

The relation (IO.30[) implies 

dX k i 



dXij 



fiikfiji + Qifokj (C.33) 



It is now straight forward to verify that 



d 



(X 2J ~% = ^Tr(X 2r )Tr(X 2J - 2r " 2 ) + Tr(X 2J - 2 ) (C.34) 



dXj ,• 

l r =0 
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Now, consider the Schwinger- Dyson equation 

o= [[dx]J-((x^) ije -^ 



which implies 

j-i 



r=0 

Writing this in terms of eigenvalues and summing over r gives 



i,j-=\,v£j 1 i=l i=l 

Now, working in terms of the eigenvalue variables we have 

JV 

= J [dx] £ (jxf- l e-^l^ (C.38) 

which becomes 

2 - 2x 2J (C.39) 



JV JV 
2J 2 2 



Comparing (1C.37j) and (1C.39I) we have 



which is solved by 



JV JV 
2 2 



2 



(C.35) 



Tr(X 2r )Tr(X 2J - 2r - 2 ) + Tr(X 2J ~ 2 )) = 2(Tr(X 2J )) (C.36) 



dl ° gJ ± + - (C.40) 



J(^) = J] a* J] (x 2 - x 2 ) = A(x) 2 (C.41) 

\k=l i<j=l J 

This is the correct answer. 

We now consider a model of two real Sp(N) matrices X and Y. Introduce the complex 
matrix Z = i (X + iY) and its Hermitian conjugate . Denote the eigenvalues of Z by 
Zi and of Z^ by Zj. Again we are interested in computing correlation functions of multitrace 
operators that have only Zs or Z^s in each trace. Given the relation between flC.6j) and 
(IC.8j) . we conjecture that the Jacobian is given by 

J(z,z) = A(z 2 )A(z 2 ) (C.42) 



35 



To prove that this guess is correct, we will again compute two point functions of our operators 
Xr(Z) and show that they agree perfectly with our result (17.211) computed directly in the 
matrix model. 

Correlation functions in the eigenvalue basis are computed using 

/N 
[dzdz}J(z, z)e~ E ^ 2 ^ • ■ ■ (C.43) 

Again fix the normalization by requiring that (1) = 1. Note that 



AT JV 
2 2 



i<j=l k=l 



^-4^+1^+1 . . . Z N ~ 2 ^ +1 (C 44) 



2 



It is now straight forward to verify that 

A/" = — ^ (C.45) 

(f)!n, 2 =i(^-2z + l)! 

Next, using the identity 

A( Z 2 )x^(Z) = c lfa "^^ + ^ +1 ^ + ^- %+ ^..z2* +W ~** +1 (C46) 

2 

we easily find 

N_ 

(X° R (Z)X°S(^)) = Uhp±"^*±W = SRS f RteveQ (C.47) 
in perfect agreement with ( I7.21j) . 
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